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Abstract: S5 0716+714 is a well-studied BL Lac object in the sky. Verifying the existence of correla- 
tions among the flux variations in different bands serves as an important tool to investigate the emission 
processes. To examine the possible existence of a lag between variations in different optical bands on 
this source, we employ a discrete correlation function (DCF) analysis on the light curves. In order to 
obtain statistically meaningful values for the cross-correlation time lags and their related uncertainties, 
we perform Monte Carlo simulations called "flux redistribution/random subset selection" (FR/RSS). 
Our analysis confirms that the variations in different optical light curves are strongly correlated. The 
time lags show a hint of the variations in high frequency band leading those in low frequency band of 
the order of a few minutes. 

Keywords: BL Lacs: general — BL Lacs: individual (S5 0716+714) 



1 Introduction 

Blazars include BL Lacertae (BL Lac) objects and fiat 
spectrum radio quasars (FSRQs). They are the most 
extreme class of active galactic nuclei (AGNs) and ex- 
hibit strong variability at all wavelengths of the whole 
electron-magnetic (EM) spectrum, strong polarization 
from radio to optical wavelengths, and are usually core 
dominated radio structures. These extreme properties 
are generally interpreted by many authors as a conse- 
quence of non-thermal emission from a relativistic jet 
orien t ed close to the line of sight (|Blandford and Konigll 
Il979l : lUrrv and Padovaiiil 1 19951 1. Blazars vary on the 
diverse time scales (Gupta et al., 2008 and references 
therein). Variability has been one of the most powerful 
tools in revealing the nature of blazars. Understand- 
ing variation behaviour is one of the major issues of 
AGNs studies. The variability lags between different 
energy bands provide very important constrains for the 
interpretation of the emission components. In recent 
years, correlations of the variabilit y in different energy 
regions have been w i dely studied ([Raiteri et al.]|200ll. 
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Ghisellini et al.iri997l:IVillata et al. 2000; Raiteri et all 



Plan et al.ll2005l : lBach et al.f 2005. 2006; N esci et al 
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2006 1: 'Montagni ct al. 2006: Fosc hini et alJ 

Wu et alj|2007l : IWagncr ct al. 1996; Stalin ct al] 



2008!: 'Dai et al."2006'; 'Mar shall et al.ll2004lArevalo et al. | 
2008; Chattcricc ct al. 2008l: IVillata et al.l2009l : lBonning e~ 
2009 ). 

The radio source S5 0716+714 is one of the bright- 
est and best-stu died BL Lacertae objects. It was dis- 
covered in 1979 (|Kiihr et al.lll98lD , and was classified 
as a BL Lac object because of its featureless spectrum 
and its strong optical polarization teiermann et al.lll98ll ). 
By optical imaging of the underlying gal axy, its red- 
shift of z = 0.31+0.08 was derived recently l|Nilsson et al.l 
120081 ). It is a highly variable BL Lac object in the 



The correlations and time lags between dif- 
ferent energy bands in this source have been stud- 
ied by some authors. A correlation was claimed for 
a selected range of data at 6 and 3 cm wavelengths , 
and optical wavelengths bv IQuirrenbach et al ] l|l99iri . 
Bu t 2 cm data from t he sa me epoch published later 
bv IQuirrenbach et all l|200.d) does not show evidence 
for correlated variability, casting doubt on the reliabil- 
ity of the cl aimed correlation b etween radio and op- 
tical bands (|Bignall et al.ll2003l ). A good correlation 
was notice d between the ligh t curves in the different 
passbands (|Sagar et al ] 1 19991 ). An upper limit to the 
possible delay between B and I varia t ions ( 10 min- 
utes) was determined by IVillata et all (|200G|) With 
the data taken on Januarv 8. 1995. lOian et al.l (l200Gl ) 
found an upper limit of the time lag of 0. 0041 days be- 
tween variations in the V and / bands. IRaiteri et al.l 
(|2003l ) found that the variations between different ra- 
alilo ba.nds are very correlated, and the flux variations 
at the lower radio frequencies are delayed with respect 
to those at the higher frequencies. They also found 
weak c orrelations b e tween the optical and radio emis- 
sions. IStalin et all l|2006l ) analyzed two night data, 
found that V and R are correlated with a small time 
lag (about 6 and 13 minutes, respectively, for April 7 
and 14, 1996) a nd the varia t ion at V leading that at R 
on both nights. IChen et al.1 (I2008D studied the gamma- 
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optical correlation. Their result suggested a possible 
delay in the gamma-ray flux variations with respect 
to optical (R band) va riations of the order of 1 day. 
iFuhrmann et al.) (|2008t ) confirmed the existence of a 
significant correlation across all their observed radio- 
bands. The time delay between the two most separated 
bands {Xsmm and Aeomm) is about 2.5 days. It must 
be noted that the time delay of 2.5 days is determined 
from an almost monotonic increase in flux density ob- 
served over a time range of 11 days. 

The determination of time lag can be used to study 
the geometry, kinematics and physical conditions in 
the inner regions of AGNs. This is attributed to light 
travel-time effect (Villata ct al. 2000). It is very im- 
portant to search for detectable delays between optical 
bands themselves, even if these can be expected to be 
very small, if any. Our goals in this paper are to inves- 
tigate the correlations between different optical light 
curves. The paper is arranged as follows: in Sect. 2, 
the light curves are presented; then in Sect. 3, the 
correlations and the time lags between different opti- 
cal passbands are presented; after this, the discussion 
and conclusions are given in Sect. 4. 



2 Light Curves 

S5 0716-f 714 is a well-monitored object. It h as been 

observed in variou s multi-frequency campaigns fWagner et al.' 
19961 : ISagar et al . 1999; ViUata ct al. 2000, 2008; Stalin c t al.. 
200d ). In addition, GhiscUini et al.1 (|1997i ) presented 



trate the fine structures of variations. Fig. [2] shows 
the light curves during JD 2,452,958-2,452,970. It is 
clear that the source does not undergo a monotonic in- 
crease between JD 2,452,958 and 2,452,970. The data 
contains many more variations than those indicated 
by Fig. [T] The source exhibits independent variation 
almost every night. 



Table 1: Statistics of BVRI-hand Data of BL Lac 
Object S5 0716+714 



Band 


N 


Mean 


a 


A 


Interval 






(mag.) 


(mag.) 


(mag.) 


(min.) 


B 


239 


14.13 


0.38 


1.42 


11.2 


V 


287 


13.69 


0.37 


1.37 


11.1 


R 


301 


13.28 


0.36 


1.37 


11.1 


I 


292 


12.75 


0.35 


1.36 


11.1 



a: 



the results of their optical ob servations b e tween Nov. 
15, 1994 and Apr. 30, 1995. iQian et al.1 (|2002l ) pub- 



lished Sra/band light curves from 1994 to 2000. iRaiteri et al] 
(|2003l ) pr esented the larges t optical database from 1994 
to 2001. IXie et al] (|2004l 'l also reporte d their moni- 
toring result on this source. [Montagn i et al.l (|2006| ) 
introduced a large set of obse rvational data contain- 
ing 10,675 photometric points. iGu et all (i2006) exten- 
sively monitored this source in 2003 and 2004. fzhang et al] 
(|2008l ) presented the optical photometries {EVRI} of 
the source from 2001 February to 2006 April. i 

To search for the lags between different optical 
bands, one needs well-sampled and high-quality datasets. 
Examining the datasets mention ed above, we fin d the 
optical photometry presented bv lGu et al.l (|2006l ) pro- 
vided a possibility to determine the time lag. The 
source was observed quasi-simultaneously in four pass- 
bands and the overall observations span about 140 
days. Each optical light curve in the BVRI bands 
is well-sampled. Furthermore, the source exhibited 
strong variability during the observations. 

The detailed statistics of these data are listed in 
Table [T] The first column represents the bands, the 
second represents the numbers, the third represents 
the mean values, the fourth represents the standard 
deviations, the fifth represents the largest variations 
and the last represents the median interval of the data 
points. The light curves of different passbands are dis- 
played in Fig. [T] They are very similar but different 
amplitudes. In the light curves, there are some signif- 
icant substructures. Because of the large time scale, 
most independent variations cannot be seen. To illus- 
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Figure 1: Light 
BVRI bands. 



curves of S5 0716+714 in the 



3 Correlation Analysis 

3.1 Discrete Correlation Function 
Method 

T he discrete correlation funct ion (DCF) was introduced 
bv lEdelson and KrolikI ||1988D . It is a useful method of 
measuring correlation, and it does not require inter- 
polating in the temporal domain. This method can 
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Figure 2: Light curves of S5 0716+714 in the 
BVRI bands during JD 2,452,958 - 2,452,970. 



not only provide the correlation of two series of un- 
evenly sampled variability data with the time lag, but 
also give the evidence of periodicity that lies in a sin- 
gle temporal data set. Its other advantages are that 
it uses all the data points available and calculates a 
meaningful error estimate. Our first step is to calcu- 
late the set of un-binned discrete correlations (UDCF) 
between each data point in the two data streams. It is 
defined as follows: 



UDCF.j = 



(oi — a) X (hj — b) 
\/<^l X 



(1) 



where Oi and hj are points of the data sets a and 
6, a and h are the means of the data sets a and b, 
(7a and at, are the standard deviations of each data 
set. Each of UDCF is associated with the pair-wise 
lag Afij = tj — ti. Then we average over the M pairs 
for which r — At/2 < Atij < t + At/2, and obtain 
DCF: 



DCF{T)^j^j:UDCF.Ar), 



(2) 



where M is the number of pairs in the bin. When 
a — b, the autocorrelation DCF is produced, and when 
a /: b, the cross-correlation DCF is measured. In most 
case, the evident peak in the cross-correlation function 
means a strong correlation between two data series, 
and the peak in the autocorrelation DCF implies a 
strong period in the data set. The standard error for 
each bin is defined as 
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Figure 3: Left: DCFs between B and /, V and /, 
R and /, respectively. Right: Normalized CCPDs 
relative to the central peak obtained by running 
8000 FR/RSS Monte Carlo simulations. Dashed 
vertical lines are drawn to guide the eye. 



a{T) = 



M ■ 



-{T.[UDCF,j - DCF{t)Y} 



2-.1/2 



(3) 



In order to obtain statistically meaningful values 
for the cross-correlation time lags and their related 
uncertainties, it is common usage to calculate the cen- 
troid Tc of the DCF, given by 



(4) 



where sums run over the points which have a DCF 
value close to the peak one {DCFi > O.SDCFpeak), 
then perform Monte Carlo simulations known as "flux 
redistribution /ran dom subset selection" (FR /RSS) de- 
scribe d in detail bv lPeterson et ahl (|l998l ) and lRaiteri et all 
(|2003l ). Random subsets of the two datasets to be cor- 
related are selected, redundant points are discarded, 
and random gaussian deviates constrained by the flux 
errors are added to the fluxes; Thus, the influence of 
both uneven sampling and flux density errors is taken 
into account. In each simulation, the two subsets are 
then cross-correlated and the centroid Tc of the DCF 
peak is determined. After a large number of simula- 
tions (generally 500~2000), the cross-correlation peak 
(actually, the centroid) distribution (CCPD) is obtained. 
As the measures of the time lag and its uncertainties, 
Tmedian and zbArgg cau be computed directly from the 
CCPD, where rfcAre s corresponds to la errors for a 
normal distribution l|Peterson et al.iri998l ). 

3.2 Result 

According to the above method, we calculate the dis- 
crete correlation function (DCF) between different op- 
tical bands to search for correlations and possible time 
lags. Between the B and I band light curves, the me- 
dian interval time is 6.3 minutes. The DCF result 
computed with a bin size of 10 minutes is shown in 
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Figure 4: Left: DCFs between B and R, V and R, 
B and V, respectively. Right: Normalized CCPDs 
relative to the central peak obtained by running 
8000 FR/RSS Monte Carlo simulations. Dashed 
vertical lines are drawn to guide the eye. 



Figure 5: Left: Auto-DCFs of B-B,V-V, R-R 
and/ — /, respectively. Right: Normalized CCPDs 
relative to the central peak obtained by running 
8000 FR/RSS Monte Carlo simulations. Dashed 
vertical lines are drawn to guide the eye. 



the top left panel of Fig O The curve of DCF has an 
obvious maximum at the position of 10 minutes. The 
maximum value of DCF is 1.1 ± 0.1. The centroid 
Tc corresponding this peak is 0.1 minutes. To give a 
statistic reliable result, we perform Monte Carlo sim- 
ulations. Because the possible delay is of the same 
magnitude as the DCF bin size and of the light curve 
time resolution, the effect of bin size cannot be ignored. 
To eliminate the effect caused by the choice of the bin 
size, we calculate the DCF with different bin sizes from 
8 to 15 minutes. The lack of a better sampling does 
not allow us to improve the bin size resolution without 
increasing spurious effects too much. We perform 1000 
Monte Carlo simulations for each bin size, thus 8000 
Monte Carlo simulations are performed. The CCPD 
is obtained and plotted in the top right panel of Fig 
|3] From this, the time lag and the uncertainties of 
10.2 



+ 13.8 



minutes are derived. In addition, the mean 
DCFpeak of 1.2 ± 0.1 is achieved by 8000 Monte sim- 
ulations. This means that the variation of B band is 
strongly correlated with that of / band with leading 
about 10.2^^3 3 minutes. 

This procedure is applied to each two bands of 
BVRI. The DCFs and CCPDs are plotted in Figs. 
13] and 2] The results are shown in Table [5] The first 
column gives the band, the second one gives the me- 
dian time interval of two bands, the third one gives the 
possible delay and the last one gives the mean value of 
DCFpeakS. The results suggest the existence of signif- 
icant correlation among BVRI band light curves and 
the higher frequencies varying earlier. 



4 Discussion and Conclusions 

It is well known that the light curves between vari- 
ant wave bands may have a short time lag. According 
to the inhomogeneous jet models, time delays are ex- 
pected between the emission in different energy bands, 



Table 2: The Possible Time Lags Determined by 
DCF 

Band Intervalimin.) Lag{min.) DCFp^ak 



B - 


I 


6.3 


10 2^20.7 


1.2 ± 0.1 


V - 


I 


4.9 


r 0-14.9 

^•°+11.4 


1.1 ± 0.1 


R- 


I 


3.7 


^■^ + 12.4 


1.0 ± 0.1 


B - 


R 


5.3 


"■'^+16.9 


1.1 ± 0.1 


V - 


R 


4.3 


Q O-10.6 
'-'■"-'+12.3 


1.1 ± 0.1 


B - 


V 


5.2 


^•"+16.8 


1.1 ± 0.1 



as pla sma disturbances propagate downstream l|Georg- anopoulos ai 
Il998l l. Multi- wavelength monitoring of blazars shows 
that flares usually begin at high frequencies and then 
propagate to lower frequencies, implying that high- 
frequency synchrotron emission arises closer to the core 
than l ow-frequency syn chrotron emission does (lUlrich et al.l 
Il997l : lMarscher|[200lh . High energy electrons emit 
synchrotron radiation at high frequencies and then cool, 

emitting at progressively lower frequencies and re sult- 

ing in time lags between high and low frequencies (iBai and Led 
(2005). The small time lags in optical regimes may be 
result of very small frequency intervals, and may indi- 
cate that the photons in these wavelengths should be 
produced by the same physical process. 

The above results have been obtained with a light- 
curve time resolution of the same magnitude as the 
possible lag. To check the method and the result men- 
tioned above, we calculate the B-B band autocorrela- 
tion DCF (Fig. EJ. The resuh of Monte Carlo simu- 
lations suggests that the time delay is — Q.O^jg g min- 
utes. For the V-V , R-R and /-/ autocorrelations, rs 



are O.O; 



-11.4 
-11. ( 



0-1+12.0 and -0.0_,_{i;g minutes(see Fig. 
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Table 3: Test Results between Light Curves and 
Their Corresponding Shifted Light Curves* 



t 


lags 


lagv 


lagR 


lag I 


1 




7-"-5 


0.6-"-^ 

-f-lZ.O 




2 


o n-13.3 
'-'■'-'+12.0 


-•^•^+11.4 


'-'•^+12.5 


^•0+13.2 


3 


1 O-12.0 
■■■■"-•+14.9 


'-'•''+9.9 


^•^+10.1 


o 4-11.7 
'-'•^+11.2 


4 


4 5-15-2 

^■^+15.0 


^•'+11.4 


^•^+11.6 


^•"+12.0 


5 


^■"+14.3 


^•"+10.8 


^•^+11.6 


4 q-11-6 

^•^+10.6 


6 


c- 7-12.8 
3- ' +12.4 


t- 7-10.8 
'+10.0 


t- 0-11.8 
^•°+11.8 


t- 0-11.9 
'^•°+11.2 


7 


7 -,-14.8 
' ■-'- + 14.4 


7Q-11.7 
' ^^+12. 6 


7 q-12.7 
' •'J+12.9 


fi q-12-2 

"•^+12.9 


8 




S 4-12.2 


R 4-12.7 


T 8-12-3 
' •"+13.3 




^•4+12.4 


^•4+12.4 


9 


q 0-12.3 
^■^+14.6 


n c;-io.o 
^•"J+ii.i 


^•^+12.4 


^■^+12.1 


10 


^^•"+12.6 


in 1-^°-^ 

^'-'•^ + 10.4 


in 

l^^"-*+11.5 


in 0-9-^ 



*The mean of each column can be seen in the 
text, and all the above are in units of minutes. 



[SJ. They are all consistent with the expectation of zero 
lags. Furthermore, we shift BVRI band light curves 
by a time of t shift and generate the corresponding ar- 
tificial light curves of B shift, VsUift, Rshift and I shift- 
Using the same process, we search for the time lag be- 
tween the light curve and the shifted light curve. The 
results are listed in Table [S] The first column is for the 
shifted time , the second to the last are for the lags of 
B-B shift, V-Vshift, R- Rshift and I-I shift determined 
by the method introduced above. From the results 
listed in Table [3l one can see that the lags between 
two correlated light curves can be well determined. In 
addition to shifting the light curves to calculate the 
displaced Auto-DCF peaks, we shift the light curves 
with respect to those at each other wavelength by the 
time lags determined. The resulting DCFs are more 
consistent with zero lag. So, the time lags between 
different optical BVRI band light curves are more re- 
liable. 

The time lag of 10.2 minutes between variations of 
B and I band is c onsistent with the result found by 
IVillata et aT](|2000r ). They derived an upper limit of 10 
minutes between B and / variations. The 5.8-minute 
lag between the band V and / is in good agree ment 
with the 0.0041-day lag determined bvlOian et al.l 12000 ') 
For the V and I band variations, IStalin et al.l (|2006l ) 
reported time lags of about 6 and 13 minutes on two 
indivi dual nights. Bet ween the V and R band varia- 
tions, [meLaD 12007') obtained the lags of 2.34 ± 5.25 
minutes on JD 2,453,737 and -0.01 ± 1.88 minutes on 
JD 2,453,742. Because the time lags are very short, 
shorter than their typical sampling interval, they con- 
cluded that they did not detect an apparent time lag. 
In this analysis, it also must be pointed out that the 
time lags between the different optical bands are of 
the same orders of the light-curve time resolution, and 
the uncertainties are larger than the lags. So, the lags 



should be treated with caution, and much more sig- 
nificant results can be achieved only with much dens er 
monit oring data. Two ways were suggested by Villata et 
(|2000l ') to improve the time resolution of the light curves, 
one way is to enlarge the telescope, another is to syn- 
chronized use two or more telescopes. 

In summary, the variations between different opti- 
cal bands have been analyzed and the time lags and 
their uncertainties are convincingly determined. The 
results suggest that the variations are correlated very 
strongly. Considering the errors, every one of the val- 
ues agrees with zero-lag well. While looking at the 
whole set of lags, it is possible to say that there is a 
hint of the variations of high frequency bands lead- 
ing those of low frequency bands in the order of a few 
minutes. 

Acknowledgments 

We express our thanks to the people helping with this 
work, and acknowledge the valuable suggestions from 
the peer reviewers. This work is supported by the Re- 
search Foundation of Education Department of Anhui 
Province, China(KJ2010B159), National Natural Sci- 
ence Foundation of China (10975145) and Fuyang Nor- 
mal University, also partly by Natural Science Foun- 
dation of Yunnan Province (2007A026M) and the In- 
novation fund (U-526) of the Institute of High Energy 
Physics. 

References 

Arevalo, P., Uttley, P., Kaspi, S., et al., 2008. MNRAS 
389, 1479. 

Bach, U., Krichbaum, T.P., Ros, E., et al., 2005. A&A 
433, 815. 

Bach, U., Krichbaum, T.P., Kraus, A., et al., 2006. 
A&A 452, 83. 

Bai, J.M., Lee, M.G., 2005. ApJ 585, L113. 

Biermann, P., Duerbeck, H., Eckart, A., et al., 1981. 
ApJ 247, L53. 

Bignall, H.E., Jauncey, D.L., Lovell, J.E.J. , et al., 
2003. ApJ 585, 653. 

Blandford, R.D., Konigl, A., 1979. ApJ 232, 34. 

Bonning, E.W., Bailyn, C, Urry, CM., et al., 2009. 
ApJ 697, L81. 

Chatterjee, R., Jorstad, S.G., Marscher, A. P., et al., 
2008. ApJ 689, 79. 

Chen, A.W., D'Ammando, F., ViUata, M., Raiteri, 
CM., et al., 2008. A&A 489, L37. 

Dai, B.Z., Zhang, B.K., Zhang, L., 2006. NewA 11, 
471. 

Edelson, R.A., Krolik, J.H., 1988. ApJ 333, 646. 



6 



Publications of the Astronomical Society of Australia 



Foschini, L., Tagliaferri, G., Pian, E., et al., 2006. 
A&A 455, 871. 

E\ihrmann, L., Krichbaum, T.P., Witzel, A., et al., 
2008. A&A 490, 1019. 

Georganopoulos, M., Marscher, A.P., 1998. ApJ 506, 
Lll. 

Ghisellini, G., Villata, M., Raiteri, CM., et al., 1997. 

A&A 327, 61. 

Gu, M.F., Lee, C.U., Pak, S., et al, 2006. A&A 450, 
39. 

Gupta, A.C., Fan, J.H., Bai, J.M., Wagner, S.J., 2008, 

AJ 135, 1384. 

Heidt, J., Wagner, S.J., 1996. A&A 305, 42. 

Kiihr, H., Pauliny-Toth, I.I.K., Witzel, A., et al., 1981. 
AJ 86, 854. 

Marscher, A.P., 2001. ASPC 224, 23. 

Marshall, K., Ryle, W.T., Miller, H.R., 2008. ApJ 677, 
880. 

Montagni, F., MascUi, A., Massaro, E., et al., 2006. 

A&A, 451, 435. 

Nesci, R., Massaro, E., Rossi, C., et al., 2005. AJ 130, 
1466. 

Nilsson, K., Pursimo, T., Sillanpaa, A., Takalo, L.O., 
Lindfors, E., 2008. A&A 487, L29. 

Ostorero, L., Wagner, S.J., Gracia, J., et al., 2006. 
A&A 451, 797. 

Peterson, B.M., Wanders, I., Home, K., et al., 1998. 
PASP 110, 660. 

Pian, E., Foschini, L., Beckmann, V., et al., 2005. 
A&A 429, 427. 

Qian, B., Tao, J., Fan, J., 2000. PASJ 52, 1075. 

Qian, B., Tao, J., Fan, J., 2002. AJ 123, 678. 

Quirrenbach, A., Witzel, A., Wagner, S., ct al., 1991. 
ApJ 372, L71. 

Quirrenbach, A., Kraus, A., Witzel, A., et al., 2000. 
A&AS 141, 221. 

Raiteri, C.M., Villata, M., AUer, H.D., et al., 2001. 
A&A 377, 396. 

Raiteri, CM., Villata, M., Tosti, G., et al., 2003. A&A 
402, 151. 

Raiteri, CM., Villata, M., Larionov, V.M., et al., 2008. 
A&A 480, 339. 

Sagar, R., Gopal-Krishna, Mohan, V., Pandey, A.K., 
Bhatt, B.C., Wagner, S.J., 1999. A&AS 134, 453. 

Stalin, C.S., Gopal-Krishna, Sagar, R., Wiita, P.J., et 
al., 2006. MNRAS 366, 1337. 



Ulrich, M.H., Maraschi, L., Urry, CM., 1997. ARA&A 
35, 445. 

Urry, CM., Padovani, P., 1995. PASP 107, 803. 

Villata, M., Mattox, J.R., Massaro, E., et al., 2000. 
A&A 363, 108. 

Villata, M., Mattox, J.R., Massaro, E., Nesci, R., et 
al., 2000. A&A 363, 108. 

Villata, M., Raiteri, CM., Larionov, V.M., et al., 2008. 
A&A 481, L79. 

Villata, M., Raiteri, CM., Larionov, V.M., et al., 2009. 
A&A 501, 455. 

Wagner, S.J., Witzel, A., Heidt, J., et al., 1996. AJ 
111, 2187. 

Wu, J., Peng, B., Zhou, X., et al., 2005. AJ 129, 1818. 

Wu, J., Zhou, X., Ma, J., et al., 2007. AJ 133, 1599. 

Xie, G.Z., Zhou, S.B., Li, K.H., et al., 2004. MNRAS 
348, 831. 

Zhang, X., Zheng, Y.G., Zhang, H.J., et al., 2008. AJ 
136, 1846. 



